##Exerccios 6 - de Teste de Hiptese e Simulao

#Crie dois conjuntos de dados: 
#A. Dez observaes de uma amostra de uma distribuio normal com mdia 6 e desvio padro 3 

help(rnorm)#rnorm(n, mean = 0, sd = 1)

conja=rnorm(10, mean = 6, sd = 3)


#B. Idem para uma distribuio normal com mdia 7.5 e desvio padro 3.2 

conjb=rnorm(10, mean = 7.5, sd = 3.2)

# 1.Utilize a funo simula.r e teste a hiptese que as mdias das amostras so diferentes. No esquea de usar a funo source() para carregar a funo!

# Estatstica de interesse - diferenas

help(source)
args(source)

#args -simula=function(dados1,dados2, nsim=1000, teste="bi")

# a funo simula j faz a diferena entre o conjunto de dados

source("http://ecologia.ib.usp.br/bie5782/lib/exe/fetch.php?media=bie5782:01_curso2009:simula.r")
simula1=simula(conja,conjb,nsim=1000,teste="bi") # ok, simulado !

# diferena absoluta observada entre as  mdias das variveis =  3.4

p1=sum(simula1>3,4)/1000
p # 0.034 - so diferentes !aceito hipotese nula

# a mdia da segunda amostra  maior que a primeira ?

simula2=simula(conjb,conja,nsim=1000,teste="uni")
simula2# 12 objetos fora 
p2=12/length(simula2)
p2 #0.012- aceito a hipotese nula

# teste t
help(t.test)
teste=t.test(conjb,conja, alternative ="greater")#unicaldal, calda com valores maiores direita
teste# 0.00515 -> aceito a hip nula- valor de p menor com teste t

# 4.Voc fez uma Anlise Exploratria dos Dados antes de fazer o teste? no =/
# est bem, ento vamos fazer agora ! da prx fao antes !

summary(conja)
summary(conjb)
str(conja)
str(conjb)
head(conja)
heas(conjb)
plot(conja)
plot(conjb)
help(boxplot)#i esqueci
boxplot(conja,conjb)#  realmente bom fazer uma anlise exploratria antes !

#5.Sempre faa um diagnstico das premissas do teste! 
# teste t- igualdade de varincias e distribuio normal 
# teste simula- as mesmas que o teste t, entretanto isso no parece ser um problema
#j que estamos usando nossos prprios dados para criar o cenrio nulo

help(bartlett.test)#Performs Bartlett's test of the null that the variances in each of the groups (samples) are the same
args(bartlett.test)


##Caixeta de NOVO?!